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ABSTRACT 

We use a high-resolution simulation of a galaxy-sized dark matter halo, published simulated 
data as well as four cluster-sized haloes from Fukushige, Kawai & Makino to study the inner 
halo structure in a A cold dark matter cos mology . We find that the circular velocity curves 
are substantially better described by Sto ehr et alJ (SWTS) profiles than by Navarro, Frenk 
& White (NFW) or Moore et al. profiles. Our findings confirm that no asymptotic slope is 
reached and that the profiles are nearly universal, but not perfectly. The velocity profiles curve 
at a constant rate in log(r) over the full converged range in radii and the corresponding ex- 
trapolated density profiles reach a finite maximum density. We find, that the claim of a strong 
discrepancy between the observed inner slopes of the density profiles of low surface bright- 
ness galaxies (LSB) and their simulated counterparts on the grounds of currently available 
observations and simulations is unfounded. In addition, if the SWTS profile turns out to be 
a good description of the halo profile for the regions that cannot be probed with simulations 
of today, then even in these regions the agreeement between simulations and observations is 
very reasonable. 

Key words: methods: N-body simulations - galaxies: clusters: general - galaxies: formation 
- galaxies: haloes - cosmology: theory - dark matter. 



1 INTRODUCTION 

One of the greatest challenges to the A cold dark matter (CDM) 
model favoured today is probably the apparent discrepancy be- 
tween the very inner parts of observed dark matter (DM) density 
profiles and their counterparts in cosmological N-body simulations 
of structure formation. 

Numerical simulations seem to find nearly 'universal' density 
profiles with cusps mostly being described by double-power-law 
fits with outer slope /? — dlogp/dlogr- ~ — 3 and inner slopes 
between -1 and -1.5. These profiles generally match the simulated 
curves reasonably well for scales of r > 0.05 r2oo where the 
virial radius r2oo is the radius of a sphere around the halo cen- 
tre enclosing a dens ity 200 times larger t han the critical density 
jNavarro. Frenk & W hite 1997, NFW; T ormen. Bouchet & Whitj 
Il997h . The exact value of the inner slope, which had been con- 
stantly re-examined as computing po wer and thus mass resolution 
increased, is still a matter of debate (Fukushige & Makino 1997 
[Moore etal. 1998, 1999b; Ghigna etal. 2000; Jine & Suto 2000 
Fukushiae & Makino 200 ll: iKlvpin et al Jl 200 1 ; Jina & Suto 2002 
iFukushige & Makino 200?; ^Power et aljbooa . The best resolved 
haloes ( Fukushiae. Kawai & Makin 



12004) with up to 3 X 10^ sim- 
ulation particles within the virial radius showed slopes at the res- 
olution limits of the simulations that were significantly shallower 
than -1.5. 

Extrapolating the fits to smaller distances results in a strong 



discrepancy with the observed rotation curves of DM-dominated 
low surface brightness (LSB) galaxies. Their D M haloes are consis- 
tent with profiles having constant density cores ( Flores & Pnmac^ 
1994; Moore 1994; Burkert 1995; McGaugh & de Blotd 1199^ 
iFirmani et al. 2001). Alth ough the extent of this discrepancy 
seemed controversial (e.g. 'van den Bosc h & Swaters 200 1.). the 
most recent studies I de Blok. McGaugh & Rubinll200lt Ide BloS 
|2004^ claim that simulations can not be reconciled with observa- 
tions. Several observational effects like seeing, misalignment of the 
slit of the telescope, finite slit size, beam smearing due to the large 
beams of the radio telescopes or an offset between the dynamical 
and the optical centre certainly matter. However, they are too small 
to make the NFW profile with inner slope of -1 mimic observed 
rotation curves ide Blok 2004). 

In this work we re-investigate the inner structure of DM haloes 
simulated at very high resolution and compare the results to pub- 
lished LSB galaxy data. In the next section we describe the simula- 
tion we have carried out and the simulated data we have used. We 
then give the results of the profile-fitting and discuss the implica- 
tions. In the last section we summarise and discuss our conclusions. 



2 SIMULATIONS 

We re ran the high-resolution simulation GA3n from lStoehr et alJ 
teOO^ with GADGET-2, a substantially improved, novel ver- 
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Figure 1. Left-hand panel: circular velocity curves for GA3new, G3-256^ fHavashi et al."2003^ and P256^ iPower et al.*2003^ (thick solid, red) with over- 
plotted SWTS (thin solid), Moore (dashed) and NFW profiles (dotted). Vertical lines indicate the radius down to which the profiles are converged. Right-hand 
panel: circular velocity curves of the four cluster simulations carried out by Fukushiee. Kawai & Makinc ( 2004). For clarity all curves have been shifted down 
by multiples of a quarter decade. The a-values of the seven SWTS profiles are 0.150, 0.128, 0.139, 0.102, 0.0934, 0.149 and 0.135. 



sioii of GADGET iSpringel. Yoshida & WhitdEooH) and termed 
it GA3new. This simulation is a two-level resi mulation start- 
ing from a large cosmologic al volume t Jenkins et alJ 1200 ll : 



I Yoshida. Sheth & Diaferij200lh carried out in a flat A-dominated 
CDM cosmology with matter density Q,m = 0.3, cosmological 
constant Q.a = 0.7, expansion rate Ha = 100 h km s~^ Mpc^^ 
with h = 0.7, fluctuation amplitude as = 0.9 at z = and box 
side length L = 479 h^^ Mpc. In the first-level resimulation, the 
mass resolution within a sphere of 52 h^^ Mpc in diameter was 
increased by a factor of 41 1. 

In order to match the properties of the Milky Way, the halo 
candidate for the GA-simulation series selected from this first-level 
resimulation was carefully chosen to be isolated and to have had a 
very quiet merging hist ory, i.e. the last ma jor merging event hap- 
pening before z — 1.5 ' Stoehr et al . 2002"). In GA3new there are 
1 1 562 566 particles within r2oo = 217 h^^ kpc at z = 0. 

The new version of GADGET combines the original tree- 
code with a particle-mesh (PM) method to compute the long-range 
forces. A 512'^-grid was used for the calculation of the correspond- 
ing fast Fourier transforms. We followed Power et al. (2003) for 
the choice of the numerical parameters but used timesteps half as 
large as required and in addition very conservative force parame- 
ters (i.e. ErrTolIntAccuracy=0.0I25, ErrTolForceAcc=0.005). The 
softening length of r^ps = 0.18 h~^ kpc was kept constant in co- 
moving coordinates. Due to the use of the PM method to compute 
the long range forces, the integration was extremely accurate also 
at high redshifts where the net forces on all particles nearly cancel 
out. 

We extracted from recently published data two additional cir- 
cular velocity curves of galaxy-sized DM haloes with lower res- 
olution. They were both produced using very similar techni ques 
as tha t described above . The haloes G3-25 6^ from Havashi et alj 
j2003h and P256^ from lPower et alj <2003h have about 2.7 x lO" 



and 3.2 X 10 particles within their virial radii, respectively. More- 
over, Toshiyuki Fukushige kindly made available the circular ve- 
locity curves of the four AC DM cluster simulations L I, L2, L3 and 
L4, presented in Fukushig e. Kawai & Makiiiol i2004) . These clus- 
ters contain 2.6 x 10^, 2.6 x lO'^, 7.2 x 10** and 7.8 x 10'^ particles 
inside their virial radii a nd were produced using GRAPE clusters 
and a GRAPE tree code jKawai & Makinj2003l) . 



3 RESULTS 

Fig. Q shows the circular velocity curves of all seven simulations 
offset by multiples of a quarter decade for clarity. Vertical lines 
indicate the radii rconv down to which the profiles have converged 
following the analysis of lPower e t al ( 2003). 

For the comparison of simulated data with analytical profiles, 
we prefer circular velocity curves over density profiles for three 
reasons. Firstly, they typically span only about half a decade on the 
ordinate, which allows to identify by eye deviations from analyti- 
cal fits in the per cent range. In addition, as they are a cumulative 
quantity, they suffer much less from Poisson noise. Finally, they are 
unique measures as no biiming is necessary. 

A possible concern due to the cumulative nature of the circular 
velocity curves is that the profiles at radii somewhat larger than the 
converged radii could be affected by the not well-resolved inner 
parts of th e haloes lead ing to underestimates of the inner slopes 
( Kazantzid is et all2004h . We will address this concern below. 

The circular velocity curves in Fig. Q (thick solid lines) do 
not asymptotically approach a constant inner slope (i.e. do not 
have an "asymptotic slope") but continue to curve all the way 
down to the converged radii. Such behaviour has been pointed 
out re c ently by several authors jSu bramanian^.Cen&_Ostrikei 
l2000j: iHavashi et all l2003l : IPower et alJ 120031 iTasitsiomi et al 
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Figure 2. Left-hand panel: standard relative deviations of the NFW (asterisks), Moore (squares) and SWTS (dots) profiles for all seven haloes. The lines 
indicate the average values (NFW: dotted, Moore: dashed, SWTS solid) excluding the values of the merging system LI (halo number 4). Right-hand panel: 
density profiles of GA3new (solid), GA2n (long dashed), GAln (dotted) and GAOn (solid) together with the SWTS density profile (smooth solid) corresponding 
to the best fit of GA3new in Fig.Q(i.e. the best-fitting combination of a, r,„aa; and Vmax)- The profiles of GA2n, GAln and GAOn, i.e. the lower resolution 
runs of the GA-simulation series, are offset by one decade for clarity. The short diagonal line indicates the slope of -1. The slope f3 of the SWTS density profile 
is shown by t he lower diagonal l ine. Its values are given on the right-hand axis. The thick lines show the average slope envelope of a recent numerical study 
carried out bv lReed et alJ i20o3l drawn down to its typical converged radius. Its radii were scaled down to match those of GA3new. The dotted fines show 
best-fit polynomials to the slope of GA3new. Vertical lines from left to right show the softening length reps, the converged radii rconv of GA3new, GA2n, 
GAln and GAOn as well as the virial radius r'200 of the GA simulation series. 



12004'; 'Hoeft, Muecket & Gottloebe^ |2004 iNavarro et all |2004 
IPiemand. Moore & StadeL^2004 '). Instead of a double-power-law 
profile, we therefore fit an analytical profile with continuously 
changing slope dlog Vc/dlogr to the profile data. The simplest 
such profile is the Stoehi: .,_Whi te. Tormen & Sorinsel (SWTS) pro- 
file JStoehretalJbooi l2003l) . i.e. a parabola, where the slope 
changes at constant rate in log(r) (Fig.Q thin solid lines): 

log {Vc/Vmax) = -a [log {r/rmax)f (1) 

Details about the corresponding mass and density profiles are 
given in the Appendix. Note that the profile fits so well, that most 
of the time, the lines are covered by the simulated data over the full 
range in r from the virial radius r2oo down to the converged radius 

^conv • 

All fits are done minimizing the standard relative deviation 

. / Er=l [(^c(n) - Ve,analyUc{ri)) /Veju)]' 
O'rel = V (2) 

V n 

measured at 100 points per decade in r. For GA3new, tJrei is only 
0.4 per cent. 

For comparis on we over plot also the best-fitting NFW 
NavMTOj Frenk & White! Il997^ (dotted) and Moore profiles 
Moore et al. 1999a) (dashed). These profiles, which have two not 
three adjustable parameters, show significant systematic deviations 
from the simulated data. For all of the seven haloes, the NFW pro- 
files under-, over- and under-predict again the data at 1 percent of 
7'200, 10 percent of r2oo and r2oo, respectively. For the Moore pro- 
files the situation is inverted. 



The left-hand panel of Fig. |2| shows for all seven haloes the 
standard relative deviations of the best-fitting NFW (asterisks), 
Moore (squares) and SWTS (dots) profiles. 

Except for the most massive cluster LI, very good SWTS fits 
are possible. The cluster LI is in a merging phase at 2 = 0: the 
slope changes abruptly at 0.4 and at 5 per cent of r2oo- In the hi- 
erarchical structure formation scenario, galaxy clusters are build- 
ing up today and are thus, by construction, subject to recent merg- 
ing events. This favours galaxy-sized DM haloes over cluster-sized 
haloes for the study of their inner structure. 

The average deviation between the data and the analytical fits 
is 4.2 per cent for the NFW profiles (dotted), 8.0 per cent for the 
Moore profiles (dashed), but only 0.89 per cent for the SWTS pro- 
files (solid). We excluded here the merging cluster LI. 

As mentioned above, the unresolved very central parts of the 
DM distribution could lead to an underestimate of the circular ve- 
locity curve at small radii. In this case, the density profile would 
have a steeper slope at small radii than the density profile corre- 
sponding to the best SWTS fit (i.e. the best-fitting combination of 
a, Vmax and Vmax) obtained from the circular velocity curve. It is 
clear from the right-hand panel of Fig. |2| that this is not the case. 
At the converged radius of GA3new, at rconv ~ 0.7 h^^ kpc, the 
SWTS density profile (smooth solid line; eQuations lAlllA4l and lA5t 
falls well on top of the density profile of GA3new (solid). 

The match is good even down to radii of about 0.5 rconv 
corresponding to roughly 2 times the softening length Tepa- This 
makes sense , as fo r GA3new, Vconv determined by the criterion of 
IPower et alJ i2003l) is a conservative limit, as the integration and 
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Figure 3. Left-hand panel: inner density profile slopes of LSB galaxies taken from lde BlolJ l2004l) who used data fromlde Blok. McGau gh & RubitJ i200lh . 
[de Blok & Bosma (2002) and Swaters et al. i2003). We converted to h ^ kpc with h = 0.7. The logarithmic density profile slopes of the analytical circular 
velocity profile fits to GA3new of Fig.0are shown as dotted (NFW), dashed (Moore) and solid (SWTS) lines. The profiles were extrapolated beyond rcon v 
of GA3new (left-hand vertical line) as indicated by the aiTow. The thin solid line shows the slope of the pa profile recently proposed by Navarro et al]i20(34l) . 
Right-hand panel: two LSB galaxy rotation curves re presenting the bulk (d ots) and the tail (fiiangles) of the shape distribution. The profiles are normalised to 
the point where the logarithmic slope equals 0.3 (see lHavashi et alj l2003t) for details). Solid lines show all seven normalised SWTS profiles of Fig. Solid 
vertical lines indicate the corresponding converged radii. 



force accuracy of GADGET has improved substantially since the 
time their study was completed. 

We have offset the density profiles of the lower-resolution sim- 
ulations of the GA convergence series (GA2n, OA In and GAOn) 
by one magnitude for clarity (see Stoehr et al. (2002, 2003) for de- 
tails). The smooth solid line shows again the SWTS density profile 
of GA3new. 

It is also obvious from the right-hand panel of Fig.|2|that dif- 
ferences between the profile and the SWTS fit as small as those that 
can be spotted in the circular velocity plots are completely invisible 
in the density profile plot due to the Poisson noise in the radial bins 
and the compressed ordinate. 

The slope (3 of the SWTS density profile is nearly constantly 
changing in log(r) over more than 2 orders of magnitude as is in- 
dicated by the diagonal line with units given on the right. At the 
resolution limit of GA3new a slope of - 1 is reached. We overplotted 
the envelope of the averaged density profile slopes as determined in 
a recent study by Reed et al. ( 2005) (thick solid lines). Their radial 
coordinate was scaled in order to match that of GA3new. 

The lower dotted lines show the best-fitting polynomials to the 
slope of the density profile of GA3new evaluated in 100 logarith- 
mically spaced bins between the converged and the virial radii. The 
polynomials (here shown are the polynomials of orders 5 to 10) 
follow very well the SWTS prediction except for very large radii. 
There is substantial scatter in the functions due to the unavoidably 
noisy slope measurement. 

When extrapolating the SWTS density profile to very small 
radii, a. flat density profile would be reached at Vfiat (see the Ap- 
pendix). However, this is out of the reach of even the best resolved 
galaxy and cluster haloes ever carried out, GA3new and L1/L2. An 



isolated DM halo would need to be simulated with approximately 
200 billion particles within r2oo in order to resolve rfiat- 

It seems reasonable to assume constant density for radii 
smaller than rfiat where the density of the SWTS profile formally 
would decrease again. 

For radii larger than r2oo> an extrapolation to the radius where 
the density equals the mean density of the Universe, seems ap- 
propriate. This radius is about 1 Mpc or approximately 5 r2oo for 
GA3new. Note however, that the SWTS profile in this region is just 
an upper bound as the density of the material outside of the accre- 
tion shock is lower than if it were virialised. 

Our findings above have important consequences for the ap- 
parent discrepancy between simulated circular velocity profiles and 
observed LSB rotation curves. 

This is because this discrepancy in fact emerges not between 
the simulations and the observations themselves, but rather between 
the extrapolated analytical power-law fits and the observations as 
has been also pointed out by Havashi et al. (2003) and de Blo3 
(2004). Extrapolation is necessary as even GA3new only is con- 
verged down to about 1 h~^ kpc whereas observationally the inner- 
most slopes of the LSB rotation curves are obtained on scales of 
0.01 to 1 h"^ kpc. 

As we have shown above, the density profile continues to flat- 
ten towards smaller radii and no asymptotic slope is reached. This 
means that slopes at different radii are necessarily different. We ar- 
gue here that much of the apparent discrepancy is a result of this 
comparison at different radii. 

The left-hand panel of Fig. |3| shows the logarithmic slopes 
of the density profiles of observed LSB galaxies as compiled by 
Ide Blokl <2004 who used data from Ide Blok. McGaugh & RubirJ 
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j200lh . Ide Blok & Bosmal <2002h and ISwaters et al] jiool. We 

overplot the NFW, Moore and SWTS density profile slopes of 
GA3new. Vertical lines indicate rconv and r2oo- 

We find that the overall agreement between the extrapolated 
SWTS profile slopes and the observed values is quite good. This 
is true even though neither a possible scatter in the a value or the 
halo mass nor observational effects like slit misplacement or beam 
smearing have been taken into account. It has been shown that these 
observational uncertainties, all resulting slope in underestimates, 
probably are too small to account for the large d iscrepancy b etween 
observed slopes and NFW or Moore profiles I de Blok 2004). How- 
ever, they may be large enough to reduce the remaining discrep- 
ancy with the SWTS profiles. This discrepancy is much smaller: 
for about half of the LSB slopes, the GA3new profile is within the 
1-(T errors. 

The thin solid line in Fig. |3] show s the slope of the ex trap- 
olated Pa profile recently proposed by iNavarro et alj i2004) that 
was fitted to the density profile of GA3new (third adjustable pa- 
rameter a = 0.18). This profile (just like the SWTS profile) does 
not reach an asymptotic slope and has been shown to be nearly 
identi cal to the SWTS pro file in the converged regions of simula- 
tions iNavarro et ai]l2004) . For the the pa profile, the agreement 
with the observed slopes is less good than that of the SWTS profile 
although it is still much closer to the observed values than it is the 
case for the NFW or Moore profiles. 

Given the necessity for extrapolation as well as the observa- 
tional uncertainities, it is clear from Fig.|3|that the claim that there 
is a strong discrepancy between the observed inner slopes of LSB 
galaxies and their simulated counterpartes does not hold. However, 
we show here in addition that, if the extrapolated SWTS profile is a 
good description of the actual DM halo profiles at small radii, then 
the agreement of observations and simulations is indeed astonish- 
ingly good. 

As a word of caution we note that, although LSB galaxies 
are DM-dominated and analysis suggests that baryons do not have 
significant influence on the inner parts of the LSB DM profiles 
I de Blok 2004), a confirmation with future simulations of LSB-type 
galaxies including realistic gas physics should be done. 

Unfortuna tely, des pite significant pr ogress made by sev- 
eral aut hors jPeeblesI fl980: Hemguis l Il99(j: IS ver & White 
| l998t iN usser & Sheth 1999; Subrama nian et all I2000bll: 
iTavlor & Navarro 200L ; Dekel et al. 2QQ.3.) , it had not been 
possible yet to determine the shape of the dark matter profiles 
based on analytical arguments. 

We now proceed from the innermost slopes to the overall 
profile shapes. The rig ht-hand panel of Fig. ID shows tw o LSB 
galaxy rotation curves ^ de Blok. McGaugh & R ubin"200l') repre- 
senting the bu lk and the tail of the LSB shape distribution (see 
rHavashietaP 12003) for details). They a re normalised to the point 
at which the logarithmic slope equals 0.3 iHavashi et alj20()3l) . The 
'typical' profile, f571-8 (dots), is very well described by all seven 
normalised SWTS profiles of Fig.Q(solid). This is true even down 
to the innermost data points for which the profiles had to be extrap- 
olated beyond their converged radii (upper vertical lines). 

The profiles of the simulated haloes fail to match the rotation 
curve of f568-3 from the tail of the profile shape distribution. The 
possible explanation, that LSB galaxies with profiles similar to that 
of f568-3 might be subject to very strong tidal fields that alter their 
inner structure, does not hold: the LSB galaxy samples contain iso- 
lated objects. 



4 SUMMARY AND DISCUSSION 

We have reinvestigated the iimer structure of DM haloes simulated 
in a ACDM cosmology using the best resolved galaxy and cluster 
haloes existing to date. The studied haloes were produced by dif- 
ferent groups using three different simulation codes. We find that 
the concept of an 'inner slope' of the circular velocity curves or 
density profiles is not appropriate: the profiles continue to curve 
without reaching an asymptotic slope in log(r-). 

The simplest analytical velocity profile with such property, the 
profile with constantly changing slope, is the SWTS profile. It de- 
scribes the analysed DM haloes typically at a level of better than 
one per cent over the full range from the converged to the virial ra- 
dius. This is substantially better than the values we obtain for NFW 
or Moore profiles. 

Allowing for one more free parameter in the analytic profiles, 
the shape parameter o, is necessary in order to provide fits of high 
accuracy. It seems plausible that the small scatter in shape parame- 
ter values reflects the differences in non-sphericity of the DM distri- 
butions. We thus find that simulated dark matter profiles are nearly 
universal but not perfectly so. Our a values range from 0.093 to 
0.15 with a median value of amed=0. 135. 

It is interesting that the same simple profile accurately de- 
scribes both the profiles of isolated haloes as well as those of 
strongly perturbed substructure haloes jStoehr et al. 2002. 2 0031) . 

There has been a claim of a strong apparent discrepancy be- 
tween the innermost density profile slopes of observed LSB galax- 
ies and their simulated counterparts. We find that this claim is un- 
founded. It relys on an extrapolation of the analytical profiles be- 
yond the actual resolution limits of the numerical simulations. Do- 
ing so is particularity dangerous, as we have shown that NFW and 
Moore profiles, which were used to support the claim, show sys- 
tematic deviations from the simulated profiles already in the nu- 
merically converged regions. 

Indeed, when extrapolating the SWTS profile (which fits the 
simulated profiles substantially better than the NFW or Moore pro- 
files) instead, we find an astonishingly good agreement with the 
innermost profile slopes of observed LSB galaxies. 

We show in addition that then the shape of a 'typical' LSB 
rotation curve can be very well reproduced. At least some of the 
remaining substantial disagreement for the shapes and slopes of 
the LSB galaxies from the tail of the distribution may be explained 
with observational effects. 
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APPENDIX A: THE SWTS PROFILE 

The density profile corresponding to the SWTS circular velocity 
from equation and its logarithmic slope (5 are 



pswTs{r) = —— IQ L V'maxyj 



47rG 



PsWTsif) 



_ 4a log (^— ) 
log [-^\ + 



(Al) 



hi{10)[l-4alog(r/r™,)] 



(A2) 



for r > Tcorii where Tcot^ is the radius at which the density profile 
becomes constant. This radius can be obtained by requiring that the 
mass interior to Tcot^ (at constant density) is the mass required to 
produce the rotation curve at rcore, i-e. 



47r 3 / 

-^rcorePSWTS(rco 
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This gives 

i_ 

The maximal density of the profile is then 

o T/2 

r \ ^ ^ max 1 r\l /(2a) 
PSWTs{rcore) = — ^ ~ 

This density is adopted for r < Vcore- The corresponding Af (< r) 
and V{r) are just (r < r-core) 



(A4) 
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<) = PSWTs{rcor 
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V{r) =r \l — G pswTs{rcor 
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It is easy to cross-check, that the core radius is the radius at which 
the rotation curve has a slope of 1. The mass Af (< r) for r > Vcore 
is 



M{< r) = 



2a [log(- 



(A8) 



G G 
We note however, that the density profile is not smooth at rcore- 
The radius at which the density profile gets flat is a bit smaller than 
rcore for typical values of a: 

I flat — ' max -LU • 

The density at this point is 

/ N 1 ^max . . . 3V° + 16°/ ln(10)-l-8a/ In(lO) 

pswTsir-fiat) = — 10 



(A9) 



3+ v/9 + 16a/ln(10) 



(AlO) 



Using rfiat instead of rcore, and thus obtaining a smooth density 
profile resulting in a smooth total rotation curve is not a bad approx- 
imation. For GA3new, the relative difference in mass integrated up 
to the maximum of the rotation curve is 4.0 x 10""". 



For GA3new, rfiat, 1' cor Si T^convi '^maxt ^200, Vmax, V200 and 
a are 0.0156 h"^ kpc, 0.0216 h"^ kpc, 0.701 kpc, 54.58 
kpc, 217.27 kpc, 251.22 km s"\ 217.27 km s"^ and 0.15, 
respectively. 



Circular velocity profiles 



